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Abstract 

The magneto-optical spectra of NiMnSb were calculated in the framework of the Local Spin Den- 
sity Approximation (LSDA) combined with Dynamical Mean-Field Theory (DMFT). Comparing 
with results based on plain LSDA, an additional account of many-body correlations via DMFT 
results in a noticably improved agreement of the theoretical Kerr-rotation and ellipticity spectra 
with corresponding experimental data. 

PACS numbers: 78.20.Ls, 71.20.Be, 71.15.Rf 
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Since the discovery of the giant magneto-optical Kerr effect (MOKE) in PtMnSb 
magneto-optical "properties became an important issue for the Mn-based family of Heusler 
alloys . However, despite of the similar structure, the group of isoelec- 

tronic alloys PtMnSb, NiMnSb and PdMnSb show quite different maximum amplitudes in 
their MOKE spectra qLMI. A theoretical description of the observed difference of the MOKE 
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spectra became possible within ab-initio band-structure calculations [8|, im [111, 113 • How- 
ever, although the various calculated MOKE spectra give reasonable qualitative agreement 
with experiment, one can notice that there exist several systematic discrepancies generally 
ascribed to the use of local spin density approximation (LSDA) . In particular, for NiMnSb, 
one can see that the low-energy peak of the Kerr rotation spectra at 1.4 eV is shifted to a 
regime of 1.6-2 eV. Also there is a noticable deviation of the amplitude for the peak at 4 eV 
as well as in the intermediate energy regime. Among other reasons, the discrepancies en- 
countered in LSDA-based results could appear due to an insufficient treatment of electronic 
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correlations. For example, as it was shown for bcc Ni |13l.ll4l|. the account of local dynamical 



correlations is essentially important for a proper description of its MOKE spectra. In the 
case of NiMnSb the main contribution to the optical transitions comes from the (i-shell of 
Mn which supplies the unoccupied part of the density of states (DOS). At the same time 
(i-electrons of Mn should be treated as locally-correlated jlS^. Based on this supposition one 
can expect an improvement of MO spectra in NiMnSb by an appropriate account of local 
correlation effects for the Mn rf-shell in band-structure calculations. In the present work the 
after is implemented within the so-called Dynamical Mean Field Theory approach (DMFT) 

In our calculations the central quantity is the optical conductivity tensor aw, as all 
optical and magneto-optical properties, can be expressed through its cartesian components 
(A = {x, y, z}). In particular for the complex Kerr an gle, that combines Kerr rotation 9k 
and ellipticity Ek, the following expression can be used [17[: 

OKico) + ieK(u;) ^ T^^^^^^l . (1) 
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In order to calculate the optical conductivity we use the expression for the hermitian com- 
ponent derived for the zero-temperature case by implementing the Green's function 
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formalism in Kubo's linear response theory [igj ]: 

a^^l(uj) = --^ y"de Tr |j^53G(e + ^a;) J^'^5G'(e)|, (2) 

where V is the volume of the spatial averaging, is the current-density operator. QG{e) 
stands for the antihermitian part of the retarded one-electron Green's function. The her- 
mitian part of the conductivity tensor 17^^(0;) can be obtained via the Kramers-Kronig 
relations. 

Formulation Q is used because it allows to include straightforwardly all correlations in 
the one-particle Green's function via the Dyson equation: 

e-i/o-S(e)]G'(e)=/, (3) 

where Hq is LSDA-based one-particle Hamiltonian including the kinetic energy, electron- 
ion interaction and the Hartree potential, while the self-energy operator S describes all 
static and dynamic effects of electron-electron exchange and correlations. The most popular 
approximation nowadays for the self-energy is DMFT which introduces it as a local, energy- 
independent exchange-correlation potential Vxc{r). As the introduction of such an additional 
potential does not change the properties of Hq we will incorporate this potential to i^LSDA 
and subtract this term from the self-energy operator. This means that the self-energy E 
describes exchange and correlation effects not accounted for within LSDA. 

The most straightforward and accurate way to solve Eq. (jHj) is to use the Korringa-Kohn- 
Rostoker Green's function method (KKR-GF) [2^. However subsequent calculations of 
optical properties within the KKR-GF method become very time-consuming due to the 
energy-dependent matrix elements of the current-density operator. A possibile alternative 
is to use the so-called variational, or fixed basis-set methods. Within such an approach the 
Green's function is represented as a sum over energy independent basis functions 



with the Green's function matrix defined as 



(2|j)e-(2|ffLSDA|j)-(^|S(e)|j) . (5) 
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The dynamical correlation correction to the one-particle hamiltonian i^LSDA is represented 
via the energy-dependent self-energy operator S(r, r ', e) which in general is a non-local quan- 
tity. Here, the self-energy is calculated via the DMFT approach that accounts only for 
ocal (on-site) correlations which are described within the Anderson impurity model (AIM) 



21[. By linking the Green's function of the effective impurity to the single-site Green's 
function derived from the k-space summation DMFT provides the self-consistent method 
to determine the self-energy accounting for on-site correlations. The actual approximation 
made in DMFT is the substitution of the non-local (k-dependent) self-energy by the single 
on-site component. The latter corresponds to the limit of the infinite coordination. Fortu- 
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nately, in many cases for 3-dimensional systems this is rather good approximation 

Although the formalism to be presented below is primarily used in connection with a 

site-diagonal self-energy, it should be stressed that any more complex self-energy can be 

used as well. In particular the formalism is able to deal with a site-non-diagonal self-energy 

occur ing for example within the GW approach [2^. 

Dealing with crystals one can make use of Bloch's theorem when choosing basic functions 
This leads to the k-dependent Green's function matrix: 



Gfjie) = [(^Jjje - (vli^LSDAljJ - (^jS(e)ljJJ . (6) 
The efficiency and accuracy of the approach is determined by the choice of {i-k). One of the 



most computationally efficient variational methods is the LMTO j2J] which allows one to 
get a rather accurate description of the valence and conduction bands in the range of about 
10 eV, which is enough for the calculation of the optical spectra {hu < 6-8 eV) . As the 
method introduces site- and angular momentum -dependent basis functions 

+ 5^ CR,;m'R'*L'(^-R'), (7) 

it perfectly fits to any single-site approximation of the self-energy. The superscripts "h" and 
"t" stand for the so-called "head" and "tail" parts of the basis functions. 

As the and {i\H\j) matrix elements in Eq. (jH)) are energy-independent it is enough 
to calculate them only once for each k-point. 
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In the framework of the DMFT the self-energy operator can be expressed in the form: 
S(f,r-",e) = ^<5RR,^e(|f-R|)y5:jf-R) 

«,RR' mm' 

X ^Zm'ir,r',e) 0(|f' - R'|) - R') , (8) 

where 0(r) = 1 if r is inside the atomic sphere and zero otherwise. Due to the special 
choice of LMTO basis functions ((Tj) only the "head" component will give the significant 
contribution to the matrix elements of the self-energy operator (jH)) leading to an extremely 
simple k-independent expression: 

ixLninxL'n) ^ J d'rd\' 
X ^tir-R)i:Zn'ir,r\e) $L(f-R). (9) 

The accuracy test of this approach was done in and the error was found to be within 5% 
which is substantially less than the approximations made for the estimation of the self-energy 
itself. 

The self-energy is calculated by using as an AIM-solver the so-called Spin-Polarised T- 
matrix plus Fluctuation Exchange (SPTF) approximation j2^l3- SPTF is a perturbative 
appoach which provides an analytical technique to sum the infinite set of the Feynmann 
diagramms for the several types of interactions in the uniform electron gas. Fortunately these 
sets of diagrams very often appear to be sufficient to describe the effects caused by dynamical 



correlations in moderately-correlated shells like Sci-electrons in transition metals 



as 



well as in systems with strong correlations j28|. This makes it a very attractive alternative 



to the so-called Quantum Monte-Carlo (QMC) technique |29| that sums all possible sets of 
diagramms of perturbation theory and therefore is much more time consuming. 

The DMFT scheme is implemented within the KKR-GF method which allows to use the 
advantages of the scattering theory formulation in the Green's function construction. For 
example, it is possible to calculate the self-energy in a self-consistent manner (parallel with 
the charge density) 0]. 

Introducing the antihermitian part of the Green's function matrix © Q'Gjj = ^{Gij — 
G*^) and taking into account translational symmetry, the hermitian part of the optical 
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conductivity (j2I) is expressed as: 



a 



(1) _ 




X 



5^J-,5(k,e)J,^'(k,e + M 



(10) 



with 



(11) 



n 



Actually the possibility to split the ^/^^(k, e) matrix elements into the energy-dependent and 
-independent parts makes the calculation of optical conductivity rather fast. 

The calculational procedure is built up as follows. The energy-dependent on-site self- 
energy S(e) is obtained from the self-consistent KKR-GF scheme j22|. For a given self- 
energy the Dyson equation (jHl) is used to obtain the effective Green's function in the basis 
set of the LMTO method. The antihermitian part '^G.ij of the effective Green's function is 
used to calculate the matrix elements of the current-density operator in Eq. The latter 
allows one to evaluate the optical conductivity given by Eq. (jTU)). Spin-orbit coupling (SOC) 
which is together with exchange splitting the actual source of MOKE is taken into account 
via the second- variation technique. 

The comparsion between the MO spectra of NiMnSb calculated within LSDA, 
LSDA-I-DMFT and the experimental results is shown in Fig. (^. The obvious conclu- 
sion is that an account of local correlations is essentially important to describe correctly 
the positions as well as the magnitudes of both low- and high-energy Kerr rotation peaks 
(situated at 1.4 eV and 4 eV). 

It is sufficient to consider in details only the real part of the MOKE spectrum, i.e. the 
rotation, as the Kerr ellipticity is a related quantity. 

The Mn d-shell indeed experiences noticable dynamical correlations as depicted by the 
self-energy plot shown in Fig. Q with the amplitudes of the imaginary component up to 
4 eV. For the Ni shell it is not so important the account of correlations as it is almost fully 
occupied. The effective Coulomb interaction is parametrised by ?7 = 3 eV and J = 0.9 eV. 
Numerical tests show that approximately the same results are obtained within the range of 
U = 3±0.5 eV. 

Considering Eqs. ()10p and pi|) we can straightforwardly analyse the modifications in the 
MO spectrum. As the matrix elements (?T-k| J^lik) are not influenced by DMFT, the two 
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FIG. 1: MO spectra for NiMnSb. Upper panel: Kerr rotation angle; lower panel: Kerr elliptic- 
ity. Solid line: LSDA calculations; broken line: LSDA+DMFT calculations. The square points 
represent the experimental results given in Refs. 3|- 

possible sources of the influence are the Green's function matrix and the change in the 
occupation numbers. However, in the present calculations the number of occupied and un- 
occupied states for each spin-projection is strictly conserved. This occurs due to the special 
construction of the self-energy represented by the so-called "particle-particle" channel 
The latter includes the channels of perturbation theory starting from the second-order and 
that is why pocessing the dominant contribution in the dynamical correlations for the d- 
electron shell. This leads in particular to the conclusion that local correlations modify only 
the interband part of the optical conductivity. The intraband contribution which is deter- 
mined by the the Green's function matrix elements at the Fermi level remains unchanged. 
Thus, one can dominantly relate the changes in the Kerr effect to the modifications in the in- 
terband contribution of optical conductivity caused by the renormalisation of the one-particle 
spectrum which is shown in Fig. (jHl). The false peak inserted by DMFT at 3 eV could be at- 
tributed first of all to the lack of the vertex corrections in our calculational scheme. In order 
to describe the optical transitions correctly, the two-particle Green's function has to be used. 
That is because the expression Q formulated in terms of the one-particle Green's function 
is an approximation which needs to be fuUfilled with the appropriate vertex corrections 3^ . 
However, the account of local correlations makes the implementation of vertex corrections in 
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FIG. 2: Spin-resolved dynamical self-energy for Mn obtained from SPR-KKR calculations using 
SPTF solver. Lower panel: real component; upper panel: imaginary component. Solid and broken 
lines correspond to majority and minority spin-components respectively. 
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FIG. 3: Spin-resolved DOS of the Mn d-shell. Solid line: LSDA calculations; broken line: 
LSDA-hDMFT calculations. 

the computational scheme rather complicated and needs a separate investigation. Another 
aspect which might be interesting to investigate is the account of nonlocal correlations. In 
conclusion, it has been demonstrated that an improved description of correlation effects on 
the basis of the LSDA-I-DMFT scheme leads to a substantially improved agreement of the 
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theoretical and experimental MO spectra. 
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